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FOREWORD 

In  this  report,  several  experiments  on  the  analysis  of  shell  structures 
using  the  STAGS-FIDAG  [5 1 computer  program  combination  are  described.  STAGS 
(STructural  Analysis  of  General  Shells)  [1]  is  a general  engineering  analysis 
computer  program  particularly  useful  for  nonlinear  problems  and  FIDAG  (Finite 
Differences  on  Arbitrary  Grids)  [5!  is  a general  program  for  finite  difference 
interpolation  on  arbitrary  grids. 

The  work  reported  here  was  supported  by  the  Air  Force  Office  of  Scientific 
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Section  1 
INTRODUCTION 

The  large  scale  engineering  analysis  program  STAGS  (STructural  Analysis 
of  General  Shells)  [1]  solves  linear  and  nonlinear  shell  problems  by 
a variational  finite  difference  technique.  It  is  a widely  used  program  known 
particularly  for  its  very  efficient  operation.  In  response  to  the  needs  of 
many  users,  a limited  variable  grid  capability  has  been  incorporated  in  STAGS 
to  treat  shell  problems  involving  rather  complex  boundaries.  The  extension  of 
this  capability  for  more  general  variable  grids  in  finite  difference  analysis 
has  been  the  main  motivation  for  the  development  of  FIDAG. 

FIDAG  is  a computer  program  for  general,  two-dimensional  functional  inter- 
polation. It  is  applicable  to  direct  interpolation  problems  such  as  that  of 
producing  interpolated  elevations  on  topological  maps,  potential  levels  in 
electrical  fields  and  displacement  and  stress  levels  in  structural  bending 
tests.  Its  primary  application,  however,  is  currently  in  the  solution  of 
partial  differential  equations,  particularly  when  it  is  advantageous  to  utilize 
non-uniform  grids. 

: 

The  FIDAG  program  has  been  continuously  evolving  over  the  past  five  years 
and,  in  its  various  forms,  has  been  applied  to  the  arbitrary  grid  finite  dif- 
ference analysis  of  boundary  value  problems  using  both  direct  and  variational 
approaches  [6],  and  a hybrid  finite  difference-element  approach  [71.  Although 
these  earlier  studies  produced  good  results,  the  full  power  of  the  program  was 
not  realized  since  the  present  Hermite  interpolation  capability  was  not  imple- 
mented. In  this  report  we  present  a number  of  results  obtained  by  linking 
FIDAG  with  STAGS  for  a more  general  finite  difference  analysis  capability. 

1.1  Program  Structure 

The  advantages  of  maximum  program  modularity  ar  well  known  and  need  not 
be  expounded  here.  When  operating  with  a program  as  complex  as  FIDAG,  it  is 
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by  far  simplest  to  have  it  operate  as  a separate,  independent  program,  communl* 
eating  with  other  programs  through  a data  base.  To  encourage  this  form  of 
utilization,  a host  program  GRIP  (GRid  Processor)  for  FIDAG  which  utilizes  a 
simple,  portable  data  base  system  for  interprogram  communication  was  developed 
and  was  used  for  this  study. 

One  option  in  GRIP  is  oriented  toward  its  use  in  the  analysis  of  elliptic 
partial  differential  equations  such  as  with  the  STAGS  (STructural  Analysis  of 
General  Shells)  program.  With  this  option,  GRIP  determines  an  astute  placement 
of  control  points  and  integration  points  as  well  as  performing  the  FID AG 
analysis.  The  terms  control  point  and  integration  point  are  defined  in 
Section  1.2  below. 

Presently,  using  the  host  program  GRIP,  the  arbitrary  grid  differencing 
capability  of  FIDAG  has  been  introduced  into  STAGS  as  illustrated  in  Fig.  1.1. 


DATA  BASE 

Fig.  1.1  Operation  of  STAGS-FIDAG  System  Using  Three 

Independent  Processors  with  Data  Base  Communication 

The  grid  generator  program  GRID  pictured  is  currently  in  a rather  limited  form 
adequate  for  test  purposes. 


2 


’St  xormr'Br* 


PHjPP 


■ ■ | — ■ 1 


wmmmmmmmmmmmmmm 


There  are  a number  of  very  compelling  advantages  to  this  organization 
including: 

• individual  processors  can  be  modified  with  minimal  side  effects 
on  other  processors, 

• other  related  processors  such  as  load  distribution  and  plotting 
programs  can  be  added  with  relative  ease, 

• efficiency  in  program  loading  is  realized  through  smaller 
"absolute"  files  for  the  programs  and 

• efficiency  for  multiple  case  studies  is  realized  since  redundant 
processes,  such  as  grid  generation,  need  not  be  repeated. 

The  efficiency  gain  referenced  in  the  last  item  above  is  particularly  important 
. in  typical  structural  analysis.  The  engineer  generally  is  interested  in  vari- 
ous static,  dynamic  and  buckling  analyses  for  a given  structure.  Once  a dis- 
cretization grid  has  been  formed  by  GRID  and  processed  for  finite  difference 
analysis  by  GRIP,  most  of  the  above-mentioned  analyses  can  be  performed  using 
STAGS  alone  with  the  information  in  the  data  base. 

1.2  Some  Definitions 

Consider  a smooth,  bounded  surface  d over  which  a continuous,  differenti- 
able function  u(x,y)  is  defined.  Let  N and  C be  two  sets  of  points,  each 
covering  d in  a fairly  uniform  fashion.  The  points  N'  will  be  called  the 
node  points  and  C the  control  points . With  each  control  point  c in  C, 
FIDAG  will  associate  a set  Nc  of  node  points  in  N which  are  near  c in 
some  sense.  Then,  for  each  control  point  c,  FIDAG  produces  a linear  trans- 
formation (coefficient  matrix)  Tc  which  facilitates  approximation  of  u, 
ux*  Uy ’ ux2>  •'*  uym  at  c terms  of  values  of  u and/or  u^  and  u at 

the  nodes  in  the  neighbor  set  N£  by  polynomial  interpolation. 

In  the  present  application  with  STAGS,  only  the  node  points  N are  con- 
structed by  the  grid  generator  GRID  along  with  an  element  table  analogous  to 
that  commonly  used  in  finite  element  analysis  systems.  The  element  table  de- 
fines a partition  of  d into  quadrilaterals  (called  EIA1 s for  Elementary  In- 
tegration Areas  in  this  report)  having  one  node  at  each  corner. 
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The  program  GRIP  then  places  a user  specified  number  of  control  points 
and  integration  points  at  quadrature  positions  in  each  EIA  for  use  in  the 
numerical  integration  of  the  energy  functional  by  STAGS. 


1.3  Overview 

The  remaining  sections  of  this  report  provide  results  from  a number  of 
tests  using  the  system  described  heretofore.  Three  basic  problems  are  treated 
in  order  of  increasing  difficulty.  In  Section  2,  rectangular  plate  bending  and 
membrane  displacement  are  analyzed  using  uniform  rectangular  grids.  Bending 
of  a clamped  disc  is  treated  in  Section  3 using  a non-rectangular  grid. 

Finally,  the  combined  effects  of  membrane  and  bending  on  a spherical  cap  are 
treated  in  Section  4 under  both  clamped  and  point  support  boundary  conditions. 
General  conclusions  about  the  operation  of  the  STAGS-FIDAG  system  in  engineering 
analysis  are  presented  in  Section  5. 
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Section  2 

RECTANGULAR  PLATE  PROBLEMS 


i 


Initial  studies  with  the  STAGS-FIDAG  processor  were  carried  out  on  a 
rectangular  plate  bending  problem  using  a uniform  rectangular  grid.  This  was 
particularly  useful  for  "wringing  out"  the  system  since  the  analytic  solution 
is  known  and  a great  deal  of  comparative  results  are  available.  Earlier 
studies  of  this  problem  have  appeared  in  [4,  5,  6,  7 1 and  elsewhere. 

In  anticipation  of  the  combined  effects  of  bending  and  membrane  behavior 
in  the  analysis  of  general  shells,  a pure  membrane  problem  was  then  studied, 
still  retaining  the  uniform  grid.  The  remainder  of  this  section  is  devoted 
to  the  results  of  these  two  studies. 


2 . 1 Plate  Bending  Problem 

The  plate  bending  problem  treated  here  is  illustrated  in  Figure  2.1. 
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Fig.  2.1  Test  Problem  for  Normal  Displacement  w 
under  Uniform  Normal  Pressure 


The  parameters  used  for  this  analysis  were: 


a = b = 0.8  in. 
h = 0.01  in 
p = 1 p.s.i. 
v = 0.3 

E = 30,000,000 
r = 0.364 


Thickness 
Normal  pressure 
Poisson's  ratio 
Young's  modulus 
12  p (1  - s»2)/E  h3  . 


The  placement  of  the  control  points  and  integration  points  was  made  using 
the  Gaussian  quadrature  placement  option  of  GRIP.  The  simplest  test  was  made 
using  a 5 by  5 uniform  grid  with  square  EIA's.  Using  this  grid,  a number  of 
tests  were  made,  varying  the  numbers  of  control  points  and  integration  points 
per  EIA  and  the  order.  The  tests  are  identified  by  the  form 


C i I j 0k 

which  means  an  i by  i array  of  control  points  and  a j by  j array  of  integra- 
tion points  were  used  in  each  EIA,  and  the  order  of  finite  difference  approxi- 
mation was  k.  Thus,  for  example,  Cl  12  03  means  that  one  control  point  was 
placed  at  the  center,  (xc>  y ) say,  of  each  EIA  and  integration  points  were 

placed  at  points  (x  -q,y  -q),  (x  -q,y  +q),  (x  +q,y  -q)  and  (x  +q, 

C C C c c c 

yc+q),  where  q = 0. 2886751346 • e and  e is  the  length  of  the  EIA  edge.  Fur- 
thermore, a cubic  approximation  was  used,  which  corresponds  to  FIDAG  coefficient 
matrices  of  order  n^  = (3  + 1)  (3  + 2)/2  = 10.  The  special  notation  0k+£  means 
that  coefficient  matrices  of  order  n^+  Z were  used.  For  example,  03  + 2 
implies  a coefficient  matrix  order  n^  + 2=12. 


In  finite  difference  analysis, constraints  on  the  trial  function  space  are 
obtained  by  collocation.  Thus,  it  appears  desirable  to  increase  the  support  of 
each  local  polynomial  without  increasing  its  degree.  A possible  approach  is  to 
combine  polynomials  on  each  EIA  by  weighted  averaging.  For  example,  in  the 
case  C3  13  03  there  are  nine,  possibly  distinct  cubics  used  on  each  EIA. 

The  use  of  only  one,  derived  from  the  given  nine,  is  clearly  a more  constrained 
system  which,  however,  no  longer  retains  the  exact  collocation  property  over 
the  entire  support.  Thus,  the  use  of  averaged  polynomials  does  not  necessarily 
produce  a trial  space  which  is  a subspace  of  that  corresponding  to  the  non- 
averaged  polynomials. 


* The  phrase  "support  of  a polynomial"  is  interpreted  as  the  entire  set  of  points 
and  function/derivative  values  utilized  in  the  definition  of  the  polynomial. 


The  test  results  corresponding  to  polynomial  averaging,  using  the  weights 
generated  for  the  numerical  quadrature  integration,  are  identified  by 

A i I j Ok  , 

where  the  interpretation  corresponds  to  that  for  Ci  Ij  Ok  described  pre- 
viously. All  of  these  results  are  summarized  in  Tables  2.1 -2.4. 

2.2  Observations  on  Plate  Bending  Results 

In  Table  2.2  we  note  that  almost  all  of  the  absolute  errors  (actual  - 
calculated)  are  negative.  We  conclude  that  the  space  of  trial  functions  de- 
fined by  the  finite  difference  process  does,  in  all  cases,  contain  functions 
outside  of  the  space  of  admissible  functions  which  yield  a lower  numerical 
energy  value  than  any  of  the  admissible  functions.  In  other  words,  we  find 
that  the  model  structure  is  more  flexible  than  the  physical  structure. 

This  property  is  in  contrast  with  conforming  finite  element  analysis  for 
which  the  trial  space  is  contained  in  the  admissible  function  space.  Thus, 
the  conforming  element  model  structure  is  less  flexible  than  the  physical 
structure.  This  property  of  conforming  finite  element  analysis  holds  generally 
whereas,  unfortunately,  the  corresponding  property  of  greater  flexibility  in 
finite  differences  does  not. 

The  weighted  mean  relative  errors  in  Table  2.4  indicate  that  all  of  the 
12  freedom  formulas  (03+2  means  (3  + 1)  • (3  + 2) /2  + 2 = 12  freedoms)  yield 

about  the  same  average  error.  The  case  A 2 1 2 0 3,  which  exhibited  a 

similar  mean  error,  indirectly  utilizes  12  freedoms  by  averaging  four  10  freedom 
cubics.  The  12  freedoms  used  for  these  cases  are  values  for  the  lateral  dis- 
placement w and  its  partial  derivatives  w^  and  w^  at  each  of  the  four 
corners  of  each  square  integration  area  (element  or  EIA)  used.  The  best  re- 
sult was  obtained  in  case  Cl  13  03+2  in  which  one  extended  (12  freedom) 
cubic  is  integrated  over  each  EIA  by  fifth  order  (3  x 3)  Gaussian  quadrature. 

We  shall  see  later  that  this  behavior  does  not  extend  to  arbitrary  .uadri- 
lateral  EIA's.  The  weights  used  for  the  mean  relative  error  calculations  in 
this  report  were  all  1. 
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TAdLt  2,1 

SQUARfc  PLATfc  RtNDlNG  PR08LLR#  5X5  GRID 
CALCULATfcO  VALUfcS  (*  10**3) 


X*.2,Y*,b 

X«,4,  Y*,8 

X»  ,6#  Y«  ,8 

x«,a,  y*,8 

actual 

3,51104 

6,31980 

8,0838/ 

6,63127 

Cl 

13 

03*2 

3,06722 

6,41545 

8,21695 

8,82774 

Cl 

13 

04 

3,71649 

6,72609 

8,61009 

9,30803 

C2 

12 

03 

3,67786 

6,03494 

8,49927 

9,13136 

C2 

12 

03*2 

3,66768 

6,41699 

8,21750 

8,82829 

C2 

12 

04 

3,01378 

6,50742 

8,35407 

b, 96668 

C 3 

13 

03 

3,69384 

6,50276 

8,33047 

8,92860 

C3 

13 

03*2 

3,66936 

b, 42193 

8,22522 

6,83324 

C3 

13 

04 

3,73611 

6,71804 

8,57150 

9,23511 

A2 

12 

03 

3,66963 

b, 41896 

8,22065 

6,8J1 3/ 

A2 

12 

04 

3,68414 

6,44622 

8,24948 

6,88024 

A3 

13 

03 

3,03886 

6,52382 

8,36673 

b,  95714 

Xs,2,  Y»,6 

X3,4, Y*,6 

X*,6»Y*,6 

X»  ,8, Y s ,6 

ACTUAL 

2,46717 

4,42208 

5,63946 

6,  k?  4961 

Cl 

13 

03*2 

2,49714 

4,47767 

5,71223 

6, 12838 

Cl 

13 

04 

2,63356 

4 ,74454 

6,09515 

6,5252/ 

C2 

12 

03 

2,56871 

4,60738 

5,8/656 

6,30706 

C 2 

12 

03*2 

2,49749 

4,47823 

5,71286 

6,12903 

C2 

12 

04 

2,63714 

4,55042 

5,80104 

6,22638 

C3 

13 

03 

2,62110 

4,52915 

5,78279 

6,20908 

C3 

13 

03*2 

2.49768 

4,47852 

5,71403 

6, 13205 

C 3 

13 

04 

2 , b30  10 

4 ,69331 

5,97946 

b , 4 0 4 4 1 

A2 

12 

03 

2,49880 

4,48051 

5,71551 

b, 13181 

A2 

12 

04 

2,61972 

4,51264 

5,77352 

6, 16773 

A3 

13 

03 

2,63795 

4,55658 

5,831 13 

6,28436 

Xs,2fY*,4 

X*,4,Y*,4 

X«,6. Y=,4 

X«,6, Y«,4 

ACTUAL 

1,44201 

2,55533 

3,23327 

3,45907 

Cl 

13 

03*2 

1,46996 

2.58452 

3,26775 

3,49509 

Cl 

13 

04 

1,64298 

2,76094 

3, 51773 

3,73985 

C2 

12 

03 

1,49662 

2,64821 

3,34/15 

3,57962 

C 2 

12 

03*2 

1.40033 

2,58508 

3,26635 

3,49568 

C2 

12 

04 

1,46762 

2,61264 

3,31441 

3,54462 

C3 

13 

03 

1 ,47263 

2,60928 

3,30135 

3,53887 

C.1 

13 

03*2 

1,46932 

2,58388 

3,26786 

3,49613 

C3 

13 

04 

1,63391 

2,70115 

3,42200 

3,64793 

A2 

12 

03 

1,46222 

2.587/3 

3,27090 

J,  498 13 

A2 

12 

04 

1,4/050 

2,60448 

3,31  149 

3,54465 

A3 

13 

03 

1 ,4b563 

2,63208 

3,36436 

3,65351 

Xs,2, Y®,2 

X».  4,Y».2 

Xe,6.Y*,2 

X«,6, Y«,2 

ACTUAL 

.49121 

,85099 

1,06225 

1,13144 

Cl 

13 

03*2 

,49892 

,86055 

1 .07171 

1 .14075 

Cl 

13 

04 

,52264 

.91317 

1 ,15346 

1,23591 

C2 

12 

03 

,51031 

.87890 

1,09358 

1 , 1637b 

C 2 

12 

03*2 

,49943 

,86102 

1,07207 

1 ,14108 

C 2 

12 

04 

,49858 

,85574 

1,07585 

1,20181 

C3 

13 

03 

,50400 

,86767 

1,07603 

1 , 14757 

C 3 

13 

03*2 

, 49806 

,86073 

1,07236 

1.14142 

C 3 

13 

04 

,51658 

,68193 

1,09550 

1 .16566 

A2 

12 

03 

,50184 

,86269 

1 ,07334 

1 , 14224 

A2 

12 

04 

,50342 

,85886 

1,06018 

1,21001 

A3 

13 

03 

,48820 

,85894 

1 , 13372 

1,26152 

""  *" 
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**** 

SQCAHb  HL  A T t 

Hfc.NO  1 NG  HKObLtM,  5X5  GRID 

« * 

* * 

AbSOLUTt  tKKG&S 

((ACTuAL  - CAICULATLD)  * 

10**3) 

• * * * 

X«,2, Y«>,8 

XS»4»Y*,8 

Xs,6,  Y«,8 

X»  ,6#  Y*,8 

Cl 

13 

03*2 

- , 04558 

-.09585 

-,1330b 

-.14647 

Cl 

13 

04 

- , 20b85 

-,40649 

-.52622 

-.62676 

C2 

12 

03 

-.10621 

-.31534 

-.41540 

-.45009 

C2 

12 

03*2 

- , 04594 

-.09639 

-,  13363 

-.14702 

C 2 

12 

04 

-.10214 

-.18782 

-.27020 

-.28561 

C3 

13 

03 

-,0b  223 

-,1831b 

- . 24660 

-.24753 

C3 

13 

03*2 

-.0.47/1 

10233 

-.14135 

-.15197 

C 3 

U 

04 

-.2234/ 

-.39844 

-, 48763 

-.55384 

A 2 

12 

03 

-.04/09 

-.099J6 

-.136/8 

-,  15010 

A2 

12 

04 

- ,07250 

-.  12562 

-.16561 

-.19897 

A 3 

13 

03 

-.  12722 

-.20422 

-,2828b 

-.27587 

X«,2,  Y«,6 

X*,4,Y*,6 

x«.6,  Y*,6 

Xs,8#Y*,6 

Cl 

13 

03*2 

-.02997 

- , 05559 

-.07277 

-.07877 

Cl 

13 

04 

-.16639 

-.32246 

-.45569 

-.47566 

C2 

12 

03 

-.10154 

-.18530 

-.23910 

-,2574b 

C 2 

12 

03*2 

-.03032 

-.05615 

-.0/340 

- ,07942 

C2 

12 

04 

- , 00997 

- . 1 28 J4 

-. 16158 

17677 

C3 

13 

03 

-.05393 

-.10707 

14333 

- , 1 5947 

C 3 

13 

03  + 2 

-.03041 

- , 0 5 b 4 4 

-.0/457 

- , 08244 

C3 

13 

04 

-.16293 

-.27123 

- , 3 4000 

-.35460 

A2 

12 

03 

- , 03  1 63 

-.05843 

-.07605 

- , 08220 

A2 

12 

04 

-.05255 

-.09056 

- , l J406 

-.13812 

A 3 

13 

03 

-.07078 

-.  13450 

1910/ 

-.234/5 

X*,2, Y*,4 

X ■ . 4 , Y • , 4 

X*  ,0, Ys  ,4 

x=,b, Y«,4 

Cl 

13 

03*2 

-.01794 

-.02919 

- , 03448 

-.03602 

Cl 

13 

04 

-.10097 

-.20561 

-.28446 

-,2807b 

C 2 

12 

03 

- , 0 5 A 5 0 

- , 092  8b 

-.11388 

-.12055 

C 2 

12 

03*2 

-.01032 

-.02975 

- ,03508 

-.03661 

C2 

12 

04 

-.02551 

-.05731 

-.08114 

-,0855b 

C 3 

13 

03 

-.03052 

-.05395 

06808 

- ,07980 

C 3 

13 

03*2 

-.01/31 

-,0285b 

- , 03459 

-,0j/t5O 

C3 

13 

04 

-,09}90 

-.  14582 

18873 

- . 1886  6 

A2 

12 

03 

-.02021 

-.03240 

- , 03/63 

-.03906 

A 2 

12 

04 

- ,02849 

-.04915 

-.0/822 

- , 08558 

A 3 

13 

03 

-.02361 

-,0767b 

-.13109 

-.19444 

X*,2/Y*«2 

X*.4, Y3.2 

Xb.g,y»,2 

X*,ri, Y*,2 

Cl 

13 

03*2 

-.0077  1 

- . 00956 

- , 00946 

-.00931 

Cl 

13 

04 

-.03143 

-.06218 

- ,09121 

-.10447 

C2 

12 

03 

-.01910 

-.02791 

-.03133 

-.03232 

C 2 

12 

03*2 

-.00822 

-.01003 

-.00982 

-.00964 

C 2 

12 

04 

-.00737 

-.00475 

-.01360 

- ,07037 

C 3 

13 

03 

-.01279 

-.01668 

-.013/6 

-.01613 

C3 

13 

03*2 

- , 00/45 

-.00974 

-.01011 

- ,00998 

C 3 

13 

04 

-.02537 

-.030-94 

-,03325 

-.03422 

A2 

12 

03 

-.0-1063 

-.01170 

-.01109 

-.01080 

A 2 

12 

04 

-.01221 

-.00787 

,00207 

- , 2/857 

* * < 

A .5 

13 

* * 

03 

,00301 

- .0G795 

-.0714/ 

-.13008 

9 


Cl  u 
Cl  u 
C2  12 
C2  12 
C2  12 
13  13 
C3  13 
C3  13 
A2  12 
A2  12 
A3  13 


03*2 

04 

03 

03*2 

04 

03 

03*2 

04 

03 

04 
03 


Cl  13 
Cl  13 
C2  12 
C2  12 
C 2 12 
C 3 13 
C3  13 
C3  13 
A 2 12 
A2  12 
A3  13 


03*2 

04 

03 

03*2 

04 

03 

03*2 

04 

03 

04 
03 


Cl  13 
Cl  13 
C2  12 
C2  12 
C2  12 
C3  13 
C 3 13 
C3  13 
A2  12 
A2  12 
A3  13 


03*2 

04 

03 

03*2 

04 

03 

03*2 

04 

03 

04 
03 


Cl  13 
Cl  13 
C2  12 
C2  12 
C 2 12 
C3  13 
C3  13 
C 3 13 
A2  12 
A2  12 
A3  13 


03*2 

04 

03 

03*2 

04 

03 

U3*2 

04 

03 

04 
03 


-1.298 
-5,891 
-4, 733 
-1,348 
-2,949 
-2,341 
-1,359 
-6,364 
-1,364 
-2,465 
-3,623 


X«,2,  Ys,6 
-1,215 
-6,744 
-4 , 1 1 1> 
-1,229 
-2,836 
-2,186 
-1,233 
-b , 6lc4 
-1,282 
-2,133 
-2,669 


X«,2,  Y*,4 
-1,244 
-7,442 
-3,784 
-1  ,27k) 
-1,769 
-2,116 
-1,240 
-6,373 
-1,441 
-1,975 
- 1 ,638 


X«,2, Y»,2 
-1,574 
-6,399 
-3,889 
-1,674 
-1,501 
-2,644 
-1.517 
-5,165 
-2,164 
-2,486 
,612 


nnuc.  c , 

N01NG  H) 
KUK  IN  ' 

. j 

<081 

rnt 

t * * < 

.CM,  5X5  GWJD 
CAlCULATfcD  VALOtS 

X*,4, Y*,8 

X3,6,Y3,8 

Xs,8, V*,8 

-1,517 

- 1 ,646 

-1,687 

-6,432 

-6,510 

-7,220 

-4,9< 

70 

-5,139 

-5,185 

-1 ,51 

2b 

- 1 , 653 

-1,694 

-2,9; 

r2 

-3,342 

-3,290 

-2,81 

78 

-3,051 

-2,851 

-1.6] 

19 

-1,749 

-1,751 

-6,31 

35 

-6 , 032 

-6,380 

• 1 • 5 1 

*2 

-1,692 

-1,729 

-1,91 

38 

-2,049 

-2,292 

-3,21 

32 

-3,499 

-3,1/8 

X * , 4 , Y 3 1 

.6 

Xs,b, Y«,6 

X*,8, Y», 6 

-1,21 

37 

-1,290 

-1,302 

-7,21 

72 

-8,080 

-7 ,863 

-4,11 

70 

-4,240 

-4,256 

-1  ,2i 

70 

- 1 ,302 

-1,313 

-2,91 

t>2 

-2,865 

-2,922 

10 


i 

TAOLt  2,4 

1# 

snoANt 

PLAffc  BLNOlUi  PKOOLfcM# 

bXb  OHIO 

* 

ML T BOOS 

SUNTLO  BY 

ML AN  KLLAflVL 

LPKOK  (D/D) 

ft  ft 

1 

C 2 

13 

03*2 

1,475 

2 

C2 

12 

03*2 

i ,4b4 

3 

A2 

12 

03 

1 , b3to 

4 

L3 

13 

03*2 

l,b4i 

5 

A2 

12 

04 

2,225 

b 

03 

13 

03 

2,724 

7 

C2 

12 

04 

3 , Dob 

0 

A3 

13 

03 

3,502 

9 

C2 

12 

03 

4 , / 08 

ID 

C3 

13 

04 

6, 109 

• 1 

2 2 

C2 

13 
* * 

04 

7,229 

ft  ft 

The  remaining  results  all  involved  more  than  the  12  freedoms  discussed 
above  and  all  gave  noticeably  poorer  results.  In  general,  it  does  not  appear 
that  significant  improvement  is  achieved  by  averaging  polynomials. 

2.3  Membrane  Problem 

The  cantilever  plate  illustrated  in  Fig.  2.3  was  analyzed  in  order  to 
obtain  purely  membrane  results.  Both  3 by  5 and  5 by  9 uniform  grids  were 
used  to  obtain  the  results  given  in  Table  2.5.  Being  a lower  order  equation 


— £ = 6 in. 
Free 


A 


E = 10,000,000 
v = • 3 
t = .1 


Free 


Symmetry 

Boundary  Cond.-Unit  vertical  displacement  at  free  end 


T 

| h = 1 in. 

Jl 


Fig.  2.3  Cantilever  Plate 


system,  only  quadratic  (02)  interpolation  was  used  in  this  test.  We  note  that 
the  averaged  polynomial  results  A 2 1 2 0 2 for  the  3x5  were  better  than  the 
non-averaged.  This  result  appears  to  occur  in  rectangular  grids  when  the 
neighbor  patterns  are  not  symmetric  about  the  EIA.  Figure  2.4  shows  the  6- 
point  neighbor  pattern  used  for  a typical  control  point  of  this  problem. 

Theoretically,  the  total  applied  force  required  to  produce  a uniform  unit 
displacement  at  the  free  end  is  given  by 

? £ ( (2  jg/h)2  + 2 (1+ v)  ) - Eht, 

which  yields  P=  8636  for  this  problem.  Thus,  we  notice  from  Table  2.5  that 
C 2 1 2 0 2 converges  from  the  flexible  side  and  A 2 1 2 0 2 from  the  stiff 
side.  In  the  remainder  of  this  report,  we  shall  see  that  the  cases  Ci  Ij  0k 
are  the  more  interesting  and  useful. 


TABLE  2.5 

Displacement  Results  along  Symmetry  Line  of  Cantilever  Plate  Problem 
and  Total  Applied  Force 


3x5  Grid 

5 x 

9 Grid 

C 2 12  0 2 

A 2 12  02 

C 2 12  0 1 

A 2 12  02 

.0758 

.0871 

.0900 

.0920 

.3065 

.3159 

.3227 

.3211 

.6488 

.6357 

.6396 

.6384 

1.0 

1.0 

1.0 

1.0 

815.0. 

9008 

8l4.42 

ms. 

Fig.  2.4  Typical  6-Point  Neighbor  Pattern 
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Section  3 
DISC  PROBIEM 


The  disc  problem  provides  an  appropriate  case  for  testing  a more  complex 
grid  structure  without  the  added  complexity  of  coupled  bending  and  membrane 
effects.  In  this  section  we  provide  results  for  disc  bending  tests  using  two 
grids.  The  conclusions  are  slightly  different  from  those  for  the  uniform  rect- 
angular grid  of  the  previous  section. 

3.1  Problem  Definition 

The  specific  problem  used  is  illustrated  in  Fig.  3.1.  The  analytic  solu- 
tion for  the  normal  displacement  in  this  problem  is  given  by 

2 2 2 2 

v(x,y)  = c (RZ-xZ-yZ)Z  , 

where 

c = 3q  (1  - v2)/a6Eh3) 

= 6.25  xlO'4 

and  R = 2 . 


Fig.  3.1  Disc  Problem 


\ 

i 

1 

I 

i 

s 

t 
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Two  grids,  designated  G and  B for  "good"  and  "bad"  were  used  for  the 
numerical  solution  to  this  problem.  Grid  B,  Fig.  3.3,  is  obtained  from  a 
square  grid  by  the  expression 


rsj 


s!M®.p  , 

S — — — 

II  Ell  2 


where  p = (x,  y)  represents  any  point  on  the  square  grid,  s is  a scale 

fSJ 

factor,  and  p'  = (x',y')  represents  the  corresponding  point  on  the  disc.  The 
norms  are  given  by 

HrII.  - max  ( |x  | , |y  |)  , 

!!  P II  2 “ V X2  + y*  . 


Grid  G,  Fig.  3.2,  was  obtained  from  B by  adjusting  the  internal  diagonal 
points  so  that  the  quadrilateral  adjacent  to  the  origin  is  square. 


Bending  results  were  obtained  for  this  problem  using  varying  numbers  of 
control  points  and  integration  points  in  each  EIA  as  described  in  Section  2.1. 
These  results  appear  in  Tables  3.1  - 3.8,  using  the  notation  introduced  in 
Section  2.1  to  identify  the  cases. 


3.2  Observations  on  Disc  Problem 

From  the  negative  signs  of  the  errors  in  Tables  3.2  and  3.6  we  note  that 
the  finite  difference  model  is  more  flexible  than  the  physical  model.  This  is 
consistent  with  the  results  for  the  rectangular  plate  discussed  in  Section  2.2. 

Of  course,  we  should  be  pleased  if  one  modus  operand i for  this  analysis 
system  could  be  shown  to  perform  well  for  any  grid.  Unfortunately,  we  have  not 
found  a theoretical  basis  for  establishing  one  and  so  we  must  look  to  experi- 
ence. If  we  let  the  mean  relative  error  be  our  guide,  we  obtain  suggestive  re- 
sults by  comparing  Tables  2.4,  3.4,  and  3.8. 

The  best  overall  performance  was  obtained  from  C 2 1 2 0 3 + 2 . For  a 
reasonably  good  grid,  Cl  13  03+2  also  performed  very  well  but  suffered 
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*car 


TAttLfc  3,1 

CLAHPtO  DISC  btNDlNG  H908LLM,  5X5  GKIO  C 
CALCULATtO  VALUtS  (*  10**2) 


Xv  1 # Y ■ 4 

X» 

2 , Y * 4 

X* 

3,Yv  4 

Xv 

4 , Y » 4 

ACTUAL 

• 1 VI  41 

,19141 

,19141 

,07368 

Cl  13  03*2 

,21012 

,21106 

,19926 

,06647 

C2  12  03*2 

,21334 

,2041b 

,20615 

,07402 

C2  12  03 

,2000/ 

,20111 

,20632 

,06926 

C3  13  03 

.2022/ 

,19713 

,19613 

,0/164 

A2  12  03 

,20/12 

,26230 

,23385 

,08323 

Xv  1 , Y ■ 3 

X* 

2 , Y ■ 3 

Xv 

3,  Y*  3 

Xv 

4 , Y • 3 

actual 

,56250 

,56250 

,40414 

,19141 

Cl  13  03*2 

,00445 

,59954 

,42652 

,19926 

C2  12  03*2 

,59624 

,60201 

,42186 

,20371 

C2  12  03 

,09652 

,60097 

,42078 

,20639 

C3  13  03 

,56966 

, 56693 

,42165 

,19661 

A2  12  03 

, /392b 

,75013 

,01613 

,23419 

Xa  1 # Y * 2 

X* 

2 » Y ■ 2 

Xv 

3,  Y v 2 

Xv 

4 , Y ■ 2 

ACTUAL 

,67891 

,76562 

,06250 

,19141 

Cl  1J  03*2 

,9297/ 

,61026 

,59954 

,21106 

C2  12  03*2 

,92560 

,00301 

,60229 

,20368 

C2  12  03 

,92014 

,60606 

,59960 

,22114 

L3  13  03 

,91741 

, 78927 

,58826 

,19/52 

A2  12  03 

1,12564 

,98672 

,77992 

,26456 

Xv  1 # Y ■ 1 

X* 

2 # Y » 1 

Xv 

3 , Y ■ 1 

Xv 

4,Yv  l 

ACTUAL 

1,00000 

,87891 

,06250 

,19141 

Cl  13  03*2 

1,05346 

.92977 

,60445 

,21 512 

C2  12  03*2 

1,04612 

,92663 

,60056 

,2115/ 

C2  12  03 

1,042/0 

,93094 

,09810 

,20494 

C3  13  03 

1,02702 

,91628 

,59014 

,22259 

A2  12  03 

1,26945 

.11967 

,77636 

,27198 

► TABLE  3,2  * 

' CLAPPED  DISC  SENDING  PROBLEM#  5X5  GRID  G * 

• absolute  errors  ((actual  - calculated)  * ib**2)  * 


X«  1,Y«  4 

X ■ 2,Y«  4 

X«  3, Y*  4 

X«  4 » Y * 4 

Cl 

13 

03*2 

-.02371 

-.01965 

-.00785 

.00721 

C2 

12 

03*2 

-.02193 

-.012/4 

-.01474 

-.00034 

C2 

12 

03 

-.01416 

-.009/0 

-.01491 

.00442 

C3 

13 

03 

- . 0 1066 

-.00572 

-.004/2 

.00184 

A2 

12 

03 

-.07671 

-.07089 

-.04244 

-.00955 

X»  1 * Y ■ 3 

X«  2.  Y*  3 

X>  3,Y>  3 

Xi  4.Y*  3 

Cl 

13 

03*2 

-.04195 

-.03704 

-.02238 

-.00/85 

C2 

12 

03*2 

-.03374 

-.03951 

-.017/2 

-.01230 

C2 

12 

03 

-.03402 

-.03847 

-.01664 

-.01498 

C 3 

13 

03 

-.02/36 

- ,02443 

-.01751 

-.00520 

A2 

12 

03 

-.17675 

-.18763 

-.11199 

-.04278 

X«  1 f V • 2 

X«  2.Y*  2 

X « 3 # Y * 2 

X«  4 . Y * 2 

Cl 

13 

03*2 

-.05086 

- , 04464 

-.03/04 

-.01965 

C2 

12 

03*2 

-.04675 

-.03739 

-.03979 

-.0124/ 

C2 

12 

03 

-.04123 

-.04044 

-.03/10 

-.009/3 

C3 

13 

03 

-.03850 

-.02365 

- , 02576 

-.00611 

A2 

12 

03 

-.24673 

-.22110 

-.21742 

-,0731b 

X ■ 1 » Y » 1 

X*  2, Y ■ 1 

X*  3,Y»  1 

X«  4 , Y * 1 

Cl 

13 

03*2 

-.05346 

-.05086 

-.04195 

-.02371 

C2 

12 

03*2 

-.04612 

-.04792 

-.03806 

-.02016 

C2 

12 

03 

-.042/5 

- .05203 

-.03560 

-.01353 

C 3 

13 

03 

-.02752 

-.03/37 

-.02764 

-.01118 

A2 

12 

03 

-,2694b 

,75924 

-.21386 

-.08057 

L8 


* * 

* * * * 

TA6LE  3,3 

**** 

ClAMRtU  DISC  bENOlNO  PRG81.tR,  bXb  GRIO 

G 

!*•< 

PERCENT  Kfc LATIVE 

ERROR  IN  THE 

CALCULATED 

VALUES 

• *** 

X«  1 , Y*  4 

X«  2,Y«  4 

X«  3#Ye  4 

X«  4 , Y ■ 4 

* Cl 

U 

0342 

-12,3b? 

-10,266 

-4,101 

9,789 

* C2 

12 

03  + 2 

-11,467 

-6 , bbb 

-7,701 

- , 4b7 

* C2 

12 

03 

-7,398 

-5,068 

-7,790 

6,003 

* C 3 

13 

03 

-6,674 

-2,988 

-2,466 

2,b01 

* A2 

12 

03 

-39,bb4 

-37,036 

-22,172 

- 1 2 , 9b7 

X a 1 # Y * 3 

X«  2 , Y ■ 3 

xs  3, Y>  3 

Xa  4 , Y * 3 

* Cl 

13 

03  + 2 

- / , 4 b8 

-b , b85 

-b , b38 

-4,101 

* C2 

12 

03*2 

-b , 998 

-7 ,024 

-4,38b 

-6,426 

* C2 

12 

03 

-6 , 048 

-6.839 

-4,117 

-7,826 

* C3 

13 

03 

*4 , 864 

-4,343 

-4,333 

-2,717 

* A2 

12 

03 

-31,422 

-33 , 3bb 

-27,711 

-22 , 3b0 

Xa  1,Y«  2 

X*  2 # Y « 2 

X*  3 , Y » 2 

X«  4 , Y » 2 

* Cl 

13 

03  + 2 

-6,787 

-b , 83 1 

-b,b8b 

-10,266 

* C2 

12 

03*2 

-3,319 

-4,884 

-7,074 

-o.blb 

* C2 

12 

03 

-4,691 

-b , 282 

-6,b96 

-b , 083 

* CJ 

13 

03 

-4,380 

-3,089 

-4,680 

-3,192 

* A2 

12 

03 

-28,072 

-28,879 

-38 , 6b2 

-38,216 

X » 1,Y«  1 

X«  2 , Y ■ 1 

X«  3 » Y ■ 1 

Xs  4 , Y a l 

* Cl 

13 

0342 

-b , 34b 

-b. 787 

-7,468 

-12,387 

* C2 

12 

03*2 

-4,012 

-b.452 

-6 , 766 

-10.632 

* C2 

12 

03 

-4,27b 

-b, 920 

-6,329 

-/ , 069 

* (.3 

13 

03 

-2,7b2 

-4.2b2 

-4,914 

-6 , 84  1 

* A2 

12 
* ** 

0J 

-2b, 94b 

6b, 38b 

-30,020 

-42,093 

T AbLt  3,4 

» * 

* clahpld 

DISC  LENDING  PHCHLhP, 

bXb  OHIO  G 

* MfcTPOOS 

SOKTfcD  8 Y 

MEAN  RELATIVE  ERROR  (0/0) 

fc  * 

* 1 

C3 

13 

03 

3,996 

a 2 

C2 

12 

03 

b,  498 

• 3 

C2 

12 

03^2 

6,714 

• 4 

Cl 

13 

0342 

6,278 

* b 

A2 

12 

03 

44,723 

k # 

TABLE  3,5 

* * * 

* *** 

CLAMHLO  DISC 

jfcNDINL  PWOHLEM# 

5X5  GHID 

b 

CALCULATED  VALUES  (*  10**2) 

1***4 

* * * ♦ 

X»  i,Y»  4 

X*  2 , Y ■ 4 X« 

3 » Y • 4 

Xa 

4 , Y ■ 4 

ACTUAL 

,19141 

.19141 

,19141 

,19141 

Cl  13  03*2 

,228/6 

,2287b 

,21732 

,20792 

C 2 12  03 

,20671 

,20048 

, 19994 

.19921 

C 2 12  03*2 

,2140S 

,20425 

,20629 

,19883 

C3  13  03 

, 20b  4 5 

,20144 

,20287 

,20261 

A2  12  03 

,25946 

,26263 

,22873 

,21216 

X«  1,Y«  3 

X»  2,Y=  3 X » 

3 » Y a 3 

X» 

4,Y«  3 

actual 

,50250 

,56250 

,56250 

.19141 

* 

Cl  13  03*2 

,60134 

,06221 

,05773 

,21732 

C2  12  03 

,00862 

,59936 

,59647 

, 19933 

• 

C2  12  03*2 

,61798 

,61243 

.60826 

,20714 

* 

C 3 13  03 

,01531 

,61058 

,60582 

,20218 

* 

A2  12  03 

,74320 

,68303 

,07106 

,22336 

* 

X«  t,Y*  2 

Xa  2, Y » 2 Xs 

3 » Y a 2 

Xa 

4,  Y*  2 

♦ 

* 

ACTUAL 

,67891 

,87891 

,56250 

,19141 

* 

Cl  13  03*2 

1 ,03387 

1.03991 

,66221 

,22875 

C2  12  03 

,93519 

,94318 

,59708 

,19957 

* 

C2  12  03*2 

,94784 

,96929 

.61204 

.20394 

* 

C 3 13  03 

,97307 

,98106 

,61195 

,2027b 

* 

A2  12  03 

1,05857 

1,06215 

,68300 

,26451 

Xa  1 , Y ■ 1 

Xa  2 , Y = 1 Xa 

3 , Y ■ 1 

Xa 

4,  Y»  1 

* 

AC  f U A 

t ,000 00 

,87891 

, 5O250 

,19141 

* 

Cl  13  03*2 

1.17421 

1 ,03387 

,06104 

,22876 

* 

C?  12  03 

1,07442 

,92728 

,60618 

,20699 

C2  12  03*2 

1,09174 

,94/51 

,61689 

,21486 

* 

C3  13  03 

1,10015 

,97 354 

,61423 

, 2&350 

* 

A2  12  03 

1,265/0 

1,04681 

,7485b 

, 25730 

* 

************* 

*«****•*, 

* * * 

********** 

* * * * 

TAbLt  3,6 

CLAPPfcD  CISC  BfcNClNC  PKOBlfcB,  5X5  OHIO  B 
ABSOLUTt  fcKfcUKS  ((ACTUAL  - CALCULATE 0)  * 10++2) 


X»  2 , Y * 4 

-.03/34 
* , 00  y 07 

-,01284 
-.01003 
- , 0/ 1 22 


3 » Y a 4 
.02591 
,00853 
,01468 
,0114b 
,03/32 


X ■ 4,Y«  4 
-,01651 
-,00780 
-,00/42 
-,01120 
-,020/b 


X « 1 , Y * 4 
-,0373b 
-,01/30 
- , 02264 
-,01404 
- , 06U0b 


Cl  13  03+2 
C2  12  03 
C2  12  03+2 
C 3 13  03 
A2  12  03 


2,Y«  3 
,099/1 
,03686 
,04993 
,04808 
, 120b3 


3 » Y « 3 
,09b23 
,0339/ 
, 0 4 b 7 b 
,04^32 
, 10856 


4 , Y * 3 
, 02b9 1 
,00/92 
, 0 1 b73 
,0  10/7 
,0319b 


1#Y«  3 

,098b4 
,04632 
, 0b548 
, 0b26 1 
, lb0/0 


Cl  13  03*2 
C2  12  03 
C 2 12  03  + 2 
C 3 13  03 
A2  12  03 


X»  2 , Y ■ 2 
-,16100 
-.0642/ 
-.09038 
-,1021b 
-,18324 


X«  3 » Y a 2 
-,09971 
-,03458 
-,04954 
-,0494b 
- , 1 2ub0 


4 , Y ■ 2 
.03734 
,00616 
,01253 
,01134 
,07310 


1 , Yb  2 
,1549b 
,05628 
,0bb93 
,094  1 6 
, 1796b 


Cl  13  03+2 
C 2 12  03 
C 2 12  03*2 
C 3 13  03 
A2  12  03 


4,  Y » 1 

,03735 
,01558 
,0?34b 
, <6  120  9 
,00589 


3,r*  i 
,09654 
,04368 
,05439 
,03  1/3 
,1660b 


X ■ 2 , Y * 1 
-,  15496 
- , 04837 
-,06860 
-.09463 
-,16790 


X«  1 , Y a 1 
-.  17421 
-,0/442 
-,091/4 
-,1001b 
-.285/6 


Cl  13  03+2 
C 2 12  03 
02  12  03+2 
C 3 13  03 
A2  12  03 


TABLE  3,7 

CLAMPED  DISC  btNDlNG  PROBLEM*  5X5  GRID  d 


PERCENT 

relative 

ERROR  IN  THE 

calculated 

VALUES 

Xs  1,Y«  4 

X*  2 , Y • 4 

Xs  3 * Y * 4 

Xs  4, Y«  4 

Cl 

13 

03*2 

-19,513 

-19,508 

-13,536 

-8,625 

C2 

12 

03 

-9,038 

-4,739 

-4,456 

-4,075 

02 

12 

03*2 

-11,828 

-6,710 

-7,7/7 

-3,676 

03 

13 

03 

-7 , 335 

-5,240 

-5,987 

-5,851 

A2 

12 

03 

-35,552 

-37,208 

-19,497 

-10,841 

x • 1,Y»  3 

xs  2 # Y ■ 3 

Xs  3 # Y s 3 

Xs  4, Y*  3 

Cl 

13 

03*2 

-17,518 

-17,726 

-16,930 

-13,536 

C 2 

12 

03 

-8,235 

-6,553 

-6,039 

-4,138 

02 

12 

03*2 

-9,863 

-8,876 

-8 , 135 

-0,218 

C 3 

13 

03 

-9,388 

-8,548 

-7,701 

-5,627 

A2 

12 

03 

-32,124 

-21,428 

-19,300 

-16.692 

x • 1,Y*  2 

Xs  2, ys  2 

Xs  3 * Y * 2 

Xs  4, Ys  2 

Cl 

13 

03*2 

-17,631 

-18,318 

-17,726 

-19,508 

C2 

12 

03 

- 6 , 4 U 3 

-7,312 

-6,148 

-4.263 

C2 

12 

03*2 

-7,843 

-10,283 

-8 , b07 

-6,546 

C 3 

13 

03 

-10 .713 

-11,622 

-8,791 

-5,924 

A2 

12 

03 

-20,441 

-20,849 

-21.422 

-38,190 

Xs  1 , Y ■ 1 

Xs  2,Y*  1 

Xs  3,Y«  1 

xs  4, Y=  1 

Cl 

13 

03*2 

- 1 7 ,421 

-17,631 

-17,518 

-19,513 

C2 

12 

03  . 

-7,442 

-5,503 

-7,766 

-6,140 

C2 

12 

03*2 

-9,174 

-7,805 

-9,669 

-12.251 

C3 

13 

03 

- 1 8 , 0 1 5 

-10,767 

-9,196 

-6,316 

A2 

12 

03 

-28,576 

-19,103 

-33,077 

-34,423 

TAbLt  3,8 

CLAMPED  DISC  BENDING  PROBLEM*  5X5  GRID  d 
METHODS  SORTED  BY  MEAN  RELATIVE  ERROR  C0/0) 

********************************************* 


C2  12  03 

02  12  03*2 

03  13  03 
01  13  03*2 
A2  12  03 


6,798 
a, 920 
9,978 
1 7 , 6138 
24,363 


greatly  (Table  3.8)  for  the  poor  grid.  Probably  the  ideal  form  would  be 
C2  13  03+2,  however,  due  to  complexities  in  the  implementation  it  has  not 
yet  been  developed.  Very  briefly,  the  complexity  involved  is  the  suitable 
correlation  of  the  2x2  array  of  control  points  with  the  3x3  array  of  in- 
tegration points  on  each  EIA  for  numerical  integration.  This  is  a development 
to  be  made  to  STAGS. 


Section  4 

SPHERICAL  CAP  PROBLEM 


The  most  difficult  test  of  this  series  involves  a shallow  spherical  cap 
requiring  consideration  of  both  bending  and  membrane  effects.  In  this  section 
we  analyze  both  an  axisymmetric  cap  problem  and  one  that  is  not.  In  the  first 
case  we  use  a converged  solution  obtained  by  the  shell  of  revolution  program 
B0S0R[2j  as  the  "actual"  solution  for  purposes  of  error  calculation.  In  the 
second  case,  we  have  no  "actual"  solution  and  simply  present  a number  of  re- 
sults from  the  STAGS-FIDAG  system,  a new  version  of  STAGS  and  a finite  element 
program  REXBAT[8]. 

4 . 1 Problem  Definition 

The  basic  problem  is  treated  as  a quarter  cap  with  symmetry  conditions 
along  the  cuts  as  illustrated  in  Figure  4.1. 

The  grid  for  this  problem  was  obtained  by  simply  projecting  the  grid  used 
for  the  disc  up  to  the  cap,  using  spherical  coordinates  a and  p . The  a , P 
coordinates  of  any  vector  defining  a point  on  the  sphere  are  the  angles  between 
the  vector  and  the  x-axis  and  y-axis  as  illustrated  in  Figure  4.1. 


We  recognize  that  this  grid  is  not  ideal  for  either  problem.  Indeed, 
since  problem  1 is  axisymmetric,  one  needs  only  a one -dimensional  grid. 
Similarly,  for  problem  2 one  needs  to  have  grid  points  concentrated  around  the 
point  supports  in  order  to  obtain  an  accurate  representation  there.  For 
problem  1 we  did  use  the  idealized  grid  with  the  BOSOR  program  to  obtain  an 
accurate  solution.  The  two-dimensional  grid  then  provided  a means  for  con- 
veniently testing  the  behavior  of  the  STAGS-FIDAG  program  on  a combined  bending- 
membrane  problem. 

Similarly,  for  problem  2 we  were  able  to  conveniently  create  a test  en- 
vironment for  the  program  which,  even  though  it  was  not  ideal  for  the  problem, 
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provided  results  indicative  of  the  program  operation.  Figure  4.2  shows  the 
appearance  of  the  grid  on  the  cap  looking  down  at  the  pole.  The  quadrilaterals 
used  in  the  STAGS-FIDAG  analysis  have  corner  points  (i,  i + 1,  i + n+l,i+n), 
where  n is  the  number  of  points  in  the  radial  direction  and  i=  (j  -l)n  + k 
for  any  j,  k=l,2,  ...  , n - 1.  The  diagonal  lines  in  Figure  4.2  on  each  EIA 
were  introduced  by  REXBAT  for  calculating  force  distribution  as  discussed 
below. 

4.2  Comparative  Finite  Element  Analysis 

A conforming  quadrilateral  element  [3]  was-  used  by  the  REXBAT  finite 
element  analysis  program  for  this  problem.  Although  this  element  is  rather  ex- 
pensive  to  use,  it  provides  very  accurate  results  which  are  particularly  useful 
for  the  evaluation  of  the  STAGS-FIDAG  system.  At  the  present  time,  the  REXBAT 
implementation  of  this  quadrilateral  element  does  not  have  a consistent  force 
distribution  capability  and  consequently,  monotonic  convergence  from  the  "stiff" 
side  cannot  be  expected  to  strictly  hold.  The  force  distribution  is  obtained 
in  REXBAT  by  dividing  each  quadrilateral  into  two  triangles  as  shown  in 
Figure  4.2,  and  using  a consistent  triangular  element  distribution. 

4.3  Numerical  Results 

The  results  of  the  STAGS-FIDAG  and  REXBAT  tests  appear  in  Tables  4.1  -4.11. 

In  the  tables,  the  point  (x,  y)  = (1,  1)  corresponds  to  the  pole,  node  1 in 
Figure  4.2.  In  general,  a point  (x,  y)  in  a table  corresponds  to  node 

1+  (x  + n*y  - 1 - n)(n-l)/4 

in  the  corresponding  grid  having  n radial  nodes.  The  grid  in  Figure  4.2  has 
n = 9. 

As  long  as  the  number  n-1  of  radial  EIA's  is  a multiple  of  4,  the  grid 
has  non-diagonal  points  which  coincide  with  those  of  the  5x5  grid  and,  conse- 
quently, those  entries  in  the  tables  correspond.  Recall,  however,  that  the 
grid  generator  adjusts  the  diagonal  nodes  corresponding  to  x = y,  x = 2,  . . . , j 

for  the  grid  improvement  and  so  the  diagonal  entries  in  the  tables  do  not 
correspond  for  different  grids. 


TABLt  4,1 

CLAMPLU  SPMfcKlCAL  CAP  PROBLtM,  6X5  GRID 

• **• 

CALCULATLO  VALUtS  (* 

10**6) 

* * 

**** 

X«  1 # Y a 4 

Xa  2 , Y • 4 

Xa  3»  Y»  4 

Xa 

4 , Y ■ 4 

* ACTUAL 

,76680 

,76b60 

,76b80 

,31480 

* Cl  13  03*2 

,86240 

, 63360 

,78621 

,27400 

* C2  12  03*2 

,Bt>222 

,79738 

,80061 

,29921 

* RtXBAT  U 

,76148 

,7671b 

,77714 

,30380 

X « 1 / Y * 3 

X s 2,Y*  3 

Xa  3, Y=  3 

Xa 

4 » Y a 3 

* ACTUAL 

2.1493B 

2,14900 

1,69700 

,76680 

* Cl  13  03*2 

2,2bb00 

2,23/70 

1,62710 

,78660 

* C2  12  U3*2 

2,24623 

2,2b660 

1 ,69890 

,80326 

* RtXBAT  U 

2,10993 

2,20180 

1,61620 

,77714 

X«  1 » Y • 2 

Xa  2 » Y ■ 2 

Xa  3, Ye  2 

Xa 

4 » Y ■ 2 

* ACTUAL 

3 ,26600 

2,87900 

2,14900 

,76680 

* Cl  13  03*2 

3,33733 

2,94810 

2,23760 

.83510 

* C2  12  03*2 

3.39680 

2,96760 

2 , 2b9b0 

,60413 

* RtXbAT  Q 

3,34220 

2,94420 

2,20180 

,76716 

X » l # Y * l 

Xa  2 , Y ■ 1 

Xa  3, Y«  1 

Xa 

4,  Y • 1 

* ACTUAL 

3,06300 

3,2bb00 

2,14900 

,/6660 

* Cl  13  03*2 

3,73210 

3,33700 

2.26760 

,86670 

« C2  12  03*2 

3,61600 

3,40120 

2,26/50 

,63248 

* KtXBAT  Q 

3,76460 

3,34220 

A************ 

2,18990 

* * 

,76148 

**** 

I 


kill. 
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♦ * 

TABLE  4,2 

k*** 

t*  *** 

**** 

CLAMPED  SPHERICAL  CAP  PROBLEM 

, 5X5  GRID 

• * 

PERCENT  RELA1IVE 

ERROR  IN  THE 

CALCULATED 

values 

**** 

X=  1 # Y»  4 

X«  2,Y*  4 

XB 

3,Y«  4 

Xb 

4 # Y ■ 4 

Cl  13  03*2 

-1 1 ,3(06 

-6,660 

-2,665 

12.961 

C2  12  03*2 

-11,265 

-4,124 

-4,546 

4,952 

REXBAT  0 

,564 

-,  1/6 

-1,481 

3,494 

Xe  1 , Ys  3 

X«  2,Y*  3 

Xa 

3,  Y*  3 

XB 

4,Y»  3 

Cl  13  03*2 

-4,9/9 

-4,128 

-1,685 

-2,566 

C 2 12  03*2 

-4,616 

-5,016 

- , 1 l 9 

-4,892 

REXBAT  Q 

-1,903 

-2,457 

-1,202 

-1,481 

X*  1 * Y ■ 2 

X>  2,Y«  2 

XB 

3 , Y « 2 

Xb 

4 * Y * 2 

Cl  13  03*2 

-2,519 

-2,400 

-4,123 

-9,049 

C2  12  03*2 

-4,356 

-3,077 

-5,142 

-5,005 

REXbAT  Q 

-2,679 

-2,265 

-2,457 

- , 1 76 

X*  1 » Y ■ 1 

X*  2 , Y ■ l 

XB 

3,  Y»  1 

Xb 

4 » Y « 1 

Cl  13  03*2 

-1,666 

-2,519 

-5,054 

-11,739 

C2  12  03*2 

-4.193 

-4.492 

-5,049 

-8,707 

♦ * 

REXbAT  Q 

-3,320 

-2,679 

k*** 

-1,903 
******  *, 

r * * * * 

,564 

*♦#« 

TABLE  4,3 

CLAMPED  SPhEKlCAL  CAP  PROBLEM,  5X5  GRID 
METHODS  SORTED  BY  MEAN  RtlATlVfc  ERROR  (070) 


1 REXBAT  Q 2,551 

2 Cl  13  03*2  3,657 

3 C2  12  03*2  4,461 
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* TABU  4,4  * 

* CLAMPED  SPHERICAL  CAP  PROBLtH,  9X9  GRID  * 

* CALCULATED  VALUES  (*  1U++6)  * 


I 


ACTUAL 
Cl  13  03*2 
C2  12  03*2 
RCXBAT  Q 

X • 1 » Y ■ 4 
,76560 
,62200 
,80640 
,7/863 

X»  2 , Y » 4 
,76660 
,82910 
,79620 
,76052 

Xi  3#  Y*  4 
,76660 
,62410 
,79490 
,78343 

X> 

4 » Y • 4 

,63300 

,64670 

.63470 

,53096 

X ■ 1,Y«  3 

X«  2»  Y*  3 

X*  3,Yi  3 

Xi 

4 # Y * 3 

actual 

2,14900 

2,14900 

1,86700 

,76660 

Cl  13  U 3+2 

2,26930 

2,2/860 

1,97890 

.62620 

C2  12  03+2 

2,2426(0 

2,23690 

1,936/0 

,79920 

KtXBAT  Q 

2,17660 

2,16060 

1,69820 

,7 8343 

X*  1 # Y • 2 

X ■ 2#  Y*  2 

X«  3#  Y*  2 

x« 

4 » Y • 2 

actual 

3,26600 

3,08600 

2,14900 

,76560 

Cl  13  03+2 

3,39760 

3,21660 

2,26220 

,63310 

C2  12  03+2 

3,37420 

3,19640 

2,25000 

,79590 

KtXBAT  Q 

3,26640 

3,11040 

2,16060 

,78062 

X « l , Y»  1 

X « 2,Y«  l 

X«  3#  Y*  l 

XI 

4»Y*  1 

actual 

3,66300 

3,25500 

2,14900 

,76560 

Cl  13  03+2 

3,79810 

3,39940 

2,26660 

,62960 

C2  12  03+2 

3,79960 

3,36940 

2,24690 

,b0610 

KtXBAT  Q 

3,69970 

3,28640 

2,1/650 

,7/863 

P ! 

♦ 

»»*****, 

TABLE  4,5 

k 1 * t 

t # * ft 

* 

CLAMPED  SPhfcKlCAL  CAP  PROBLEM 

, 9X9  BRIO 

♦ 

• * 

PERCENT  RELATIVE 

ERROR  IN  THE 

CALCULATED 

VALUES 

: ' 

* 

* 

X« 

1 , Y»  4 

Xi  2 , Y • 4 

Xa 

3 , Y » 4 

Xa 

4 , Y • 4 

* 

Cl  13  03*2 

-7,339 

-6,266 

-7,613 

-2,303 

[ ** 

* 

C2  12  03*2 

-b ,302 

-3,9/0 

-3,800 

-,319 

e 1 • 

♦ 

RExbAT  Q 

-1 ,675 

-1,922 

-2,302 

,303 

If  ! 

* 

X* 

1 » Y ■ 3 

X*  2 , Y ■ 3 

X* 

3 # Y ■ 3 

Xa 

4,  Y«  3 

« 

Cl  13  03*2 

•b , 1 33 

-6,031 

-4,6/0 

-7,/b/ 

I 

* 

C2  12  03*2 

-4 , 3b 1 

-4,090 

-2 ,634 

-4.361 

( 

* 

REXBAT  Q 

- 1 ,280 

-1,460 

-,  594 

-2,302 

• 

* 

X« 

1 , Y ■ 2 

X»  2, Y«  2 

X> 

3 , Y » 2 

Xa 

4 # Y « 2 

I 

♦ 

Cl  13  03*2 

-4,376 

-4,200 

-6,198 

-8,706 

t 

* 

L 2 12  03*2 

-3,602 

-3,642 

-4,700 

-3,931 

* 

KEXbAT  Q 

-1  ,(426 

-,791 

-1,460 

-1,922 

fi 

* 

x« 

1 » Y*  1 

Xa  2 , Y ■ 1 

X a 

3 » Y ■ 1 

Xa 

4 , Y » 1 

1 

• 

Cl  13  03+2 

-3,668 

-4,436 

-b  , 43b 

-8,30/ 

r • ) 

♦ 

C2  12  03+2 

-3,729 

-4,129 

-4 ,009 

-0,202 

»? 

# 

♦ 

* * 

FEXbAT  U 

-1 ,002 

-1,026 

k * * ♦ 

-1,260 

*******4 

* * * * 

-1,6/0 

* TABLE  4,6 

* CLAPPED  SPHERICAL  CAP  PROBLEM,  9X9  BRIO 

* METHODS  SDK  f 1 0 BY  MEAN  RELATIVE  ERROR  (070) 


KEXbAT  Q 
C2  12  03*2 
Cl  13  03*2 


lfM6 

3,900 

4,898 


TAdLfc  4,7 

CLAMPfcD  SPhtKICAL  CAP  PRObLfc*>  13X13  GRID 
CALCULATfcD  VALUfcS  (*  10**6) 


X s l,Va  4 

X"  2,Y«  4 

X * 3 » Y • 4 

* # * 

XI 

4,Y«  4 * 

ACTUAL 

,76580 

,76580 

,76580 

.61240  • 

Cl  13  03*2 

,81280 

,82130 

,82830 

,03900  * 

RfcXBAT  G 

, 7b  1 46 

,78230 

,78365 

,61277  * 

X • 1 , Y * 3 

X«  2, Y*  3 

Xs  3 » Y ■ 3 

X * 

4, Y«  3 * 

ACTUAL 

2,14900 

2,14900 

1,98300 

,76580  * 

Cl  13  03*2 

2,25650 

2,28310 

2,09640 

,82880  * 

RfcXBAT  Q 

2,17320 

2,17520 

1,98860 

,7836b  * 

Xs  1 , Y«  2 

X * 2,Y*  2 

X«  3, Y«  2 

X* 

4,  Y*  2 * 

ACTUAL 

3, 23500 

3,14700 

2,14900 

,76580  ♦ 

Cl  13  03*2 

3,40690 

3,30100 

2,28340 

,822/0  * 

REXBAT  0 

3,27480 

3,16010 

2, 17520 

,78230  * 

X«  1,Y«  1 

X:  2,Y*  1 

X«  3 » Y s l 

X ■ 

4 » Y * 1 * 

ACTUAL 

3,66300 

3,25500 

2, 14900 

,76580  * 

Cl  13  03*2 

3,82020 

3,40680 

2,25850 

.81610  * 

# * 

# + 4 

RfcXBAT  Q 

3, b6050 

3,27480 

2,17320 

* + + 
* * * 

,/6146  * 

TAbLt  4,8 

CLAPPfcO  SPHERICAL  CAP  PKt)BL*>,  13X13  GRID 
PfcHChNT  H t L A T I V R £KHUR  IN  Thfc  CALCULATfcD  VALUfcS 


***♦♦*♦♦** 


Cl  13  03*2 
RfcXBAT  0 

X * 

1,Y«  4 
-6,137 
-2,045 

X ■ 

2,Y«  4 
-7,247 
-2,155 

X * 

3 * Y » 4 
-8,161 
-2,331 

X« 

4 , Y ■ 4 
-4,344 

- ,060 

Cl  13  03*2 
RfcXBAT  0 

x« 

1,Y»  3 
-5,002 

-1,126 

X* 

2 » Y 8 3 
-8,240 
-1,219 

X« 

3 , Y * 3 
-5,719 
-.282 

X 8 

4#  Y8  3 
-8,227 
-2,331 

Cl  13  03*2 
HtXBAT  Q 

x> 

l r Y ■ 2 
-4,667 
-,608 

x» 

2,Ys  2 
-4,894 
-,416 

Xs 

J,Y*  2 
-0,254 
-1,219 

Xs 

4 , Y * 2 
-7,430 
-2,155 

Cl  13  03*2 
RfcXBAT  Q 

Xs 

1 # Y • l 
-4,292 
-.478 

X • 

2,Y«  1 
-4 , 664 
-,608 

Xs 

3 , Y • 1 
-5,W9b 
-1,126 

Xs 

<»fYs  1 
-6,568 
-2,045 

TABLt  4,9  * 

CLAPPfcO  SPhfcRICAL  CAP  PROBLfc*,  13X13  GRID  * 

Mt T HODS  SOKTfcO  BY  MEAN  RtLATIVfc  fcRRUR  (1570)  * 

l RfcXBAT  0 ,882  * 

2 Cl  13  03  + 2 0,134  • 

V 

> 

1 

i 

TAbLE  4,10 

a 

n 

SPHERICAL  CAP, 

4 PUINT  SOPPUHT#  9X9  GRID 

| 

* 

calculated  values  (* 

10**5) 

** 

Xa  1 , Y • 4 

X«  2,Y»  4 

Xa  3 # Y « 4 

• *4 
Xa 

4,  Y*  4 

**** 

* 

* 

Cl  13  03*2 

,46/50 

,62720 

,01900 

,04920 

* 

* 

C2  12  03*2 

,48774 

,66020 

, b5610 

,09940 

• 

i 

1 ! 
i 

REXBAT  0 

,39644 

,66120 

,72520 

,76103 

♦ 

* 

X ■ 1 * Y • 3 

X*  2 , Y ■ 3 

Xa  3,Y«  3 

Xa 

4,Y«  3 

♦ 

Cl  13  03*2 

,0/400 

,94090 

,94970 

,01670 

* 

C2  12  03*2 

,94390 

,96990 

,9/900 

,06700 

1 

♦ 

REXBAT  Q 

,7/601 

,03633 

,0409b 

,72520 

* 

Xa  1 , Y*  2 

X»  2 , Y * 2 

Xa  3,Y«  2 

Xa 

4,Y«  2 

* 

1 

Cl  13  03*2 

1,40630 

1,07090 

,93910 

,62340 

# 

I 

* 

C2  12  03*2 

1 . 1 1670 

1 ,10200 

,9/630 

,66010 

♦ 

f 

* 

REXBAT  Q 

,97630 

,96099 

,03533 

,55120 

* 

0 

1 

* 

X*  l,Ya  1 

X*  2 , Y ■ 1 

X«  3 » Y*  1 

X« 

4 , Y ■ 1 

: 

•) 

* 

Cl  13  03*2 

1 ,14410 

1,00560 

,07220 

, 40460 

* 

* 

C2  12  03*2 

1,17794 

1 ,12090 

,90000 

,49150 

* 

J 

* 

♦ 

REXbAT  G 

1,03300 

.9/630 

,77601 

,39044 

y 

I 


I 


r #******  + *6  + 1 

TABLE  4,11 

SPHERICAL  CAP,  4 POINT  SUPPORT,  13X13  GRID 
CALCULATED  VALUES  <*  10**5) 


* * 

,********#■*,*< 
Xa  1 , Y a 4 

X«  2,Y«  4 

♦ * * * 
x« 

3*  Ya  4 

♦ * * 
X a 

4,  Y»  4 

« # * * 

♦ 

* 

Cl  13  03*2 

, 5»960 

,67370 

,07460 

,91460 

* 

RLXbAT  Q 

,44070 

,60452 

,78290 

, 022t0 1 

* 

X=  1 , Y a 3 

X«  2 # Y « 3 

Xa 

3,Ya  3 

X» 

4 » Y * 3 

♦ 

* 

Cl  13  03*2 

.94120 

1,01150 

1 

,03104 

,0/300 

* 

REX0AT  Q 

,04707 

,90764 

,92276 

,70298 

* 

Xa  t,Ya  2 

Xa  2, Y«  2 

Xa 

3 » Y a 2 

X« 

4 f Y ■ 2 

Cl  13  03*2 

l , 16640 

1 ,16020 

1 

,01020 

,6/110 

* 

REXbAT  Q 

1,05500 

1,04000 

,90764 

,60452 

* 

Xa  1 , Y • 1 

X*  2 , Y a 1 

Xa 

3,Ya  1 

Xa 

4,  Y»  1 

* 

Cl  13  03*2 

1,22930 

1 ,16570 

,93940 

,50730 

* 

REXbAT  Q 

1,11640 

1,05500 

,04/07 

* * * 

,440/0 

* 

* 

* + ** 
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Based  on  experience  with  the  plate  and  disc  problems,  methods  Cl  13  03  + 2 
and  C 2 1 2 0 3+f  were  used  in  this  series  of  tests  for  the  bending.  We 
recognize  that  the  general  quartic  integration  (13)  is  beneficial  and  that 
more  than  one  control  point  per  EIA  (C  2)  is  also  beneficial.  It  appears  that 
the  ideal  combination  for  accuracy  would  be  C 2 13  03+2  which  has  not  yet 
been  implemented. 

4.3.1  Clamped  Cap  Problem 

The  results  for  the  clamped  cap,  problem  1,  appear  in  Tables  4.1  -4.9 
along  with  "actual"  results  which  were  obtained  from  BOSOR  as  discussed  earlier. 

Notice  in  Tables  4.2,  4.5  and  4.8  that  the  REXBAT  solutions  do  not  converge 
from  the  "stiff"  side,  i.e.  , the  coarse  grids  produce  a solution  which  is 
greater  than  the  actual  at  most  points,  as  was  anticipated  from  the  discussion 
in  Section  4.2.  From  these  tables  we  also  notice  that  the  STAGS-FIDAG  model  is 
predominantly  more  flexible,  which  is  consistent  with  the  results  of  the  plate 
and  disc  studies.  We  shall  observe  a reversal  of  this  trend  in  the  results  for 
problem  2 below.  We  notice  in  Tables  4.6  and  4.9  that  both  STAGS-FIDAG  results 
and  the  REXBAT  results  have  essentially  converged  with  the  9x9  grid  to  mean 
errors  of  57.  and  1%,  respectively.  Also,  we  observe  from  the  similarity  of  the 
results  for  the  two  STAGS-FIDAG  forms  that  both  are  reasonably  good  and  neither 
is  clearly  superior.  This  is  suggestive  of  a probable  benefit  from  utilizing 
the  best  of  both  forms  via  C 2 1 2 0 3+  2 as  mentioned  earlier. 

One  surprising  observation  is  the  occasional  appearance  of  a slight  in- 
crease in  mean  relative  error  in  the  STAGS-FIDAG  results  when  the  grid  is  re- 
fined. For  example,  compare  the  results  for  Cl  0 3 0 3+2  in  Tables  4.3, 

4.6  and  4.9.  All  we  can  conclude  from  these  results  is  that  convergence  is 
not  monotonic. 
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4.3.2  Point  Support  Problem 

Our  results  for  this  problem  are  less  detailed  than  previous  results 
since  the  true  solution  is  not  known.  In  Tables  4.10  and  4.11  we  note  that 
all  of  the  solutions  increase  as  the  grid  is  refined.  This  is  traditional  be- 
havior for  finite  elements  but  is  unusual  for  finite  differences.  It  is  the 
only  case  in  this  study  for  which  the  finite  difference  solutions  converge  from 
the  stiff  side. 

Since  a thorough  convergence  study  was  beyond  the  budgetary  scope  of  this 
effort,  further  analysis  was  conducted  with  STAGS  as  reported  in  Section  4.4. 

In  addition,  crude  projections  based  on  REXBAT  results  from  the  5x5,  the 
9x9  and  the  13  x 13  grid  were  made.  These  appear  in  Figure  4.3.  Both  studies 
suggest  that  the  13  x 13  STAGS-FIDAG  result  is  close  to  being  correct. 

4 . 4 STAGS  Analysis 

The  grid  described  in  Section  4.2  cannot  be  used  in  STAGS  without  FIDAG. 

The  approach  used  for  a spherical  cap  or  disc  with  STAGS  involves  a "geographic" 
coordinate  system  of  longitudes  and  latitudes.  Since  the  longitudes  coalesce 
at  the  pole,  it  is  necessary  to  perturb  the  problem  by  introducing  a small 
hole  at  the  pole  and  applying  suitable  boundary  conditions  there.  The  best 
results  obtained  to  date  have  been  with  free  edge  conditions  around  the  hole. 

The  cap  considered  in  this  report  corresponds  to  a maximum  latitude  of 
about  14.5°.  For  the  STAGS  study,  the  polar  opening  used  was  1°  latitude. 

The  errors  in  the  STAGS  C'  solution  for  the  clamped  cap  are  given  in 
Table  4.12  along  one  longitude.  The  STAGS  C'  results  for  the  four-point 
support  cap  along  a longitude  of  a support  are  interpolated  to  the  latitudes 
corresponding  to  the  results  of  Table  4.11  and  presented  with  those  results 
in  Table  4.13. 
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Radial  Nodes 


Projected  Solution  from  R'L'BAT  Results  for  4 Point 
Support  ap  at  Various  Latitudes  along  a Longitude 
of  Support 


TABLE  4.12 

Errors  in  the  STAGS  C ' Results 
for  the  Clamped  Cap  at  Various  Latitudes  (x  10^) 


Section  5 
CONCLUSIONS 


5 . 1 Cost  Factors 

The  efficiency  of  the  STAGS  program  is  one  of  its  outstanding  features  and 
is  well  known  to  its  many  users.  Thus,  discussion  here  is  directed  to  the  im- 
pact to  STAGS  of  incorporating  the  arbitrary  grid  finite  difference  capabili- 
ties of  FIDAG. 

As  noted  in  Section  1,  the  FIDAG  analysis  of  the  grid  is  carried  out 
separately  from  the  STAGS  analysis,  and  thus  can  be  pro-rated  over  studies  in- 
volving a variety  of  boundary  conditions.  From  STAGS  point  of  view,  instead  of 
generating  finite  difference  coefficients  internally,  it  simply  reads  them  from 
the  data  base.  Consequently,  for  typical  engineering  analysis,  particularly  in 
nonlinear  studies,  the  cost  difference  between  treating  n uniformly  placed 
nodes  and  n arbitrarily  placed  nodes  is  very  minor.  However,  since  the  ar- 
bitrary placement  generally  allows  adequate  coverage  of  a structure  with  sub- 
stantially fewer  nodes,  the  use  of  arbitrary  grids  will  generally  result  in  a 
considerably  lower  analysis  cost. 

For  a general  grid  on  a shell,  FIDAG  is  required  to  produce  two  coefficient 
matrices  per  control  point:  a low  order  one  for  membrane  and  a higher  order 

one  for  bending.  The  results  of  this  study  indicate  that  the  FIDAG  cost  per 
control  point  using  8x8  membrane  and  12  x 12  bending  matrices  is  about  0.08 
seconds  on  the  CDC  6600.  With  6x6  membrane  and  10  x 10  bending  matrices 
the  cost  is  about  .05  seconds  per  control  point. 

5.2  Results 

The  sequence  of  tests  presented  in  this  report  constitute  a preliminary 
evaluation  of  the  STAGS-FIDAG  system.  Bending  and  membrane  problems  have  been 
treated  separately  and  in  combination  via  shell  analysis.  The  results  have 
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provided  a basis  for  establishing  suitable  value  ranges  for  the  parameters  in- 
volved, namely:  (1)  the  number  of  control  points  per  element,  (2)  the  order 

of  the  numerical  integration  to  be  used  (number  of  integration  points  per 
element)  and  (3)  the  order  of  finite  difference  approximation  to  be  used. 
Further  refinement  of  these  ranges  will  be  accomplished  through  broader 
problem  experience. 

In  general  applications,  one  will  usually  obtain  conservative  results  with 
the  system,  i.e.,  predicted  displacements  will  tend  to  be  larger  than  the 
theoretically  exact  ones. 
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